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Abstract 

We consider nucleon-decay-like signatures of the hylogenesis, a variant of the an- 
tibaryonic dark matter model. For the interaction between visible and dark matter 
sectors through the neutron portal, we calculate the rates of dark matter scatterings 
off neutron which mimic neutron-decay processes n —> 1/7 and n —> ve + e~ with richer 
kinematics. We obtain bounds on the model parameters from nonobservation of the 
neutron decays by applying the kinematical cuts adopted in the experimental analy¬ 
ses. The bounds are generally (much) weaker than those coming from the recently 
performed study of events with a single jet of high transverse momentum and missing 
energy observed at the LHC. Then we suggest several new nucleon-decay like processes 
with two mesons in the final state and estimate (accounting for the LHC constraints) 
the lower limits on the nucleon lifetime with respect to these channels. The obtained 
values appear to be promising for probing the antibaryonic dark matter at future un¬ 
derground experiments like HyperK and DUNE. 
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1 Introduction 


Given a variety of spatial scales and cosmological epochs associated with dark matter phe¬ 
nomena, their natural explanation seems in introducing a new neutral particle, stable at 
cosmological time-scales. Many extensions of the Standard Model of particle physics (SM) 
suggest suitable dark matter candidates with masses ranging from 10” 23 eV (oscillating scalar 
held, see e.g. [18]) to 10 16 GeV (superheavy dark matter, see e.g. [13]). Dark matter particles 
must be produced in the early Universe at a stage before matter-radiation equality. Most 
mechanisms exploited for this purpose work properly (for a review see [2]) but treat the 
(order-of-magnitude) equality of dark matter and visible matter contributions to the present 
energy density of Universe, 

Pdm,o ~ Pb,o , (1) 

as an accidental coincidence. 

Yet it may be a hint towards a common origin of both cosmological problems, dark matter 
phenomena and matter-antimatter asymmetry of the Universe. There are models addressing 
this issue. In particular, an elegant approach is provided by models of antibaryonic dark 
matter, where dark matter particles carry (anti)baryonic charge. The idea is that the total 
baryonic charge of the Universe is zero, but it is redistributed between visible sector (positive 
baryonic charge) and dark sector (negative charge of the same amount). Both dark and 
visible matter emerge during the same process at some stage in the early Universe making 
a connection between the two components, so that the coincidence (1) may be understood. 

Similar to the visible sector, the dark sector is asymmetric, being populated solely with 
particles of negative baryonic charge. The models of this type are called asymmetric dark 
matter, for a review see [21]. They exhibit quite specific phenomenology. As a rule, no 
dark matter pair annihilation is expected in galaxies or inside the Sun (see, however, [15, 
3]). Instead, the antibaryonic dark matter particle may annihilate with nucleon, mimicking 
proton/neutron disappearance or decay. 

A remarkable example of the antibaryonic dark matter model is provided by hylogene- 
sis [9]. To the SM particle content at low energies the model adds a complex scalar <f> and 
Dirac spinor T, together forming dark matter components, and also two heavy fermions X a , 
a = 1,2, playing the role of messengers between the visible and dark sectors. The interaction 
terms read 

C = -^-X~ a cf • tfC 1 ±Jl d k + CaXAb C <F + h.c., (2) 

A z 2 2 

with i,j, k running over the SM three generations, d l and v? denote down-type and up-type 
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quarks, superscript C refers to charge conjugation; and ( a are dimensionless coupling 
constants, A stands for the scale of new physics which completes the model to a renormal- 
izable theory (for a particular variant of high-energy completion within a supersymmetric 
framework see [6]). 

The new fields carry baryonic charge, so that B(X a ) = 1 and P(T) = P(<f>) = —1/2. 
Coupling constants X l J k and ( a are, in general, complex numbers, providing the model with 
charge (C) and charge-parity (CP) violation required for the successful dynamical generation 
of the baryon asymmetry. The latter is produced in the early Universe via CP-violating 
decays of nonrelativistic messengers X a in a way very similar to what happens in the standard 
leptogenesis with heavy sterile neutrinos [12], Since the baryon number is conserved by 
interactions (2), in the same process the dark sector (T, <f>) becomes asymmetric, collecting 
the negative baryonic charge produced in the CP-violating decays of nonrelativistic fermions 
X a . Later in the LIniverse, baryons and antibaryons of the visible sector annihilate, leaving 
the net baryonic charge, which is accumulated at present mostly in hydrogen and helium. A 
similar process happens in the dark sector, and the antibaryonic charge of the same amount 
is distributed between fermions T and bosons <f>. This may be characterized by a ratio of 
their present number densities, 


V = 


ny,o 
r> 'V ,o 


(3) 


Proton and both dark matter particles, T and <f>, are stable, if their masses obey the 
kinematical constraints 


| My — My | < Mp + m e < My + My , 


(4) 


where M p and m e stand for proton and electron masses. Total baryon number conservation 
implies a simple relation between dark matter and visible baryon number densities 


n B 


Tty Tly 
2 


(5) 


For the present dark matter energy density, one can write 


PDM ,0 — X4y Uy + A4y Tly . 


( 6 ) 


Without any asymmetry between the two dark matter components, i.e., when 
obtain from (6) and (5) 

XIy + My 


Pdm,o 


M p 


PB ,0 


n 


1, we 


(7) 
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which for the observed property (1) settles the dark matter mass scale in the GeV-range. 
Then, for the present cosmological estimates of pb,o and Pdm ,o [20], the sum of the dark 
matter particle masses are fixed by eq. (7), while the kinematical constraint (4) confines the 
individual masses inside the interval 

1.7 GeV < M# , M* < 2.9 GeV. (8) 

With asymmetry between T and populations, rj ^ 1, the relation (7) is replaced with 

2(Afjr + rj M$) 

PDM -° ~ (1 + ,,)M„ PB '° ' (9) 

The interaction with quarks in (2) can be used to probe the model at colliders [10, 11]. 
Heavy fermions X a can be directly produced or virtually contribute to dark matter produc¬ 
tion. This model provides the following signatures for the LHC experiments (depending on 
the quark structure in (2)): (%) missing energy and either a jet with high transverse momen¬ 
tum p T [10, 11] or a heavy quark (t,b, or c) with high p T [11]; (ii) a jet (or a heavy quark) 
with high px and a peak in the invariant mass of three jets whose momenta compensate 
high px [11]. The performed analysis of LHC events with a high-p^ jet and missing energy 
has allowed us to constrain the model parameter space pushing the new physics up to TeV 
scale [11], 

Another very pronounced signature of the model [9] is an induced nucleon decay (IND) [10]. 
The dark matter particle scattering off a nucleon (through the exchange of virtual fermions 
X a ) flips its type, T o <f>, and destroys the nucleon. The kinematical constraint (4) ob¬ 
viously forbids the traceless disappearance of the nucleon, i.e., a process like $ + n —> 4/. 
Some additional particles must emerge in the final state yielding a signature of the induced 
nucleon decay. These processes involving an additional single meson in the final state have 
been analyzed [9, 10, 6] for a set of quark operators entering (2) and a number of final states. 
While the scattering mimics the nucleon decay, the kinematics of particles in the final state 
is different, which prevents us from direct use of the limits on the proton/neutron lifetimes 
to constrain the model parameter space. However, by adjusting properly the kinematical 
cuts, the corresponding analysis has been performed [6, 10]. In particular, for X a couplings 
to the uds operator in (2), the results of nucleon decay searches raise the mass of heavy 
fermion X a and the scale of new physics A up to the TeV scale [9, 10, 6]. 

In this paper we analyze several new modes of the induced nucleon decays via a neutron 
portal, represented by dud operator in (2). The paper is organized as follows. In Sec. 2 
we derive the low energy effective lagrangian describing the dark matter scattering off a 
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neutron and give the relation between the scattering cross section and the nucleon lifetime 
with respect to decay into a given final state. In Sec. 3 we consider 2 —)• 2 scattering processes 
v h(<h)+n —y < h( v h)+ 7 , which mimic neutron decay n —y 1 / 7 , and, imposing the cuts adopted in 
the experimental search for this decay mode [5, 19], we constrain the model parameter space. 
These constraints turn out to be (much) weaker than those following from the LHC [11], so 
finally we obtain a lower estimate of the neutron lifetime in this model based on the limits 
from the LHC. In a similar way, we investigate the scattering T( < h) + n —> < h(T) + e + e~ 
in Sec. 4. We study the induced nucleon decays into two light mesons (7r, K, 77 in various 
possible combinations) in Sec. 5 and (based on the LHC bounds [11]) predict the shortest 
lifetimes at the level of 10 32 — 10 33 yr expected for these modes within hylogenesis. The 
obtained numbers are quite promising and allow the processes to be tested with the next 
generation underground facilities like HyperK [16, 1] and DUNE [17]. We expect that these 
channels apart from the dominant single-meson-induced nucleon decays would be helpful 
to discriminate between different models predicting processes with baryon number violation 
and corner an interesting region in the parameter space of the hylogenesis scenario if a signal 
of nucleon-decay-type is found in future. We conclude in Sec. 6 . 


2 Low-energy effective lagrangian and nucleon lifetime 


The coupling terms in eq. (2) relevant for low-energy phenomenology of the neutron portal 
read 

c = -^x~ a d ■ PP d + + h. c .. ( 10 ) 

A- 2 2 

Hereafter we are interested in processes with typical energies much below the mass scale 
of the heavy fermions X a . The exchange of virtual X a between the visible sector and dark 
sector fields entering ( 10 ) yields the following contact interaction 

v 2 xd °- udQ 1 , , , 

_ La 1 M Xa l+Jl d . U C l+Jl d + h . c . . (11) 

A 2 2 2 v ; 

For further analysis it is convenient to introduce variables Mx and y by relations 

y = y-Arc 

M x M Xa 

so that M x (somewhat vaguely) indicates the heavy fermion scale, while dimensionless pa¬ 
rameter y reflects the coupling strength. The physical meaning of M x is the energy scale 
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below which the effective interaction (11) can be safely exploited instead of (10). Since not 
y and Mx individually but only their ratio ( 12 ) enters all the formulas below, there is an 
ambiguity in the definition of y and Mx related to the change of the variables. However, it 
has no impact on the physical observables. 

Further, the GeV scale of dark matter masses ( 8 ) and smallness of the expected velocity 
of galactic dark matter particles allow us to describe the dark matter scattering off nucleons 
in terms of baryons and mesons rather than quarks and gluons. In this approximation, the 
lagrangian (11) with replacement (12) transforms into Yukawa-type interaction 

r yP asific 1 i i viqs 

£ = ~A+ < 13 ) 

which we use below to calculate the scattering rates; n denotes the neutron held and the 
parameter f3 = 0.012 GeV 3 is related to the QCD scale [ 8 ]. 

The cross sections of dark matter scatterings off nucleon N, and cr$ N ^ ..., are 

related to the total nucleon lifetime with respect to a particular IND process tjv-*.... as follows 

Tiv-k.. = --, (14) 

VCT\^x—y... T n<j) var$x^y... 

where v is the dark matter particle velocity in the laboratory frame where nucleons are 
at rest. In fact, since the scatterings we discuss happen in s-wave, the cross sections are 
inversely proportional to v, and the lifetime (14) does not depend on its value. 

3 Scattering processes ^(&)n —> <£>(\I /)7 

We start our study with a simple 2 —» 2 scattering with dark matter particles annihilating 
a neutron into dark matter particle of another type and a photon. Let p^, p n and q be the 
4-momentum of T, neutron n, and the outgoing photon 7 , being real for 4/ n —> $ 7 , and 
hence q 2 = 0 (or virtual for T n —> $e + e“, which we consider in Sec. 4). The process is 
proceeded due to the Yukawa interaction (13) and the neutron dipole moment 

£ = na ^q l ' F 2(q 2 )nA tl , (15) 

where e M (g) is photon polarization 4-vector and for the Pauli (magnetic) form factor we 
utilize the dipole parametrization F 2 (q 2 ) = — 1 . 91/(1 + q 2 r 2 M /12) 2 with magnetic radius 
r M = 0.86 fm [ 20 ], 
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The dark matter particle scatters off the neutron by means of virtual neutron exchange. 
The matrix element of the process reads 


iey/3 F 2 (q 2 ) 
2M n A 2 M x 




1 + 75 P-M n 

2 p2_ M 2 


a^q^p^e^q)®, 


where p = q — p n is the 4-momentum of the virtual neutron, and n(p n ), + (pm ), <f> are wave 
functions of the neutron, T and <f> particles, respectively. In the laboratory frame the neutron 
is at rest, while the dark matter particle moves with small velocity v <C 1. Here and below, 
we perform the estimates to the leading order in velocity v. The squared matrix element 
averaged over spins of the two incoming fermions in the laboratory frame is 


\M\ 2 


eV/? 2 Ff(0) 

4 Ml M\ A 4 


M n Mm . 


(16) 


For the similar process <f> n —> T 7 , we find the same expression (16) up to the following 
replacement 


M n Mq, —> 2M n (M$ + M n — q 0 ), 


where the additional factor accounts for different numbers of fermions in the initial states 
averaged over spins. To the leading order in v « 1, the photon frequency is 


q 0 ~ M n + Mm — Mm . 

We can place a bound on the model parameter space from nonobservation of the decay 
n —> 1/7 [5, 19] exhibiting the same signature as the scattering process under discussion: a 
single photon in the final state. To this end, we constrain the kinematics of the photon as it 
has been adopted 3 in the original experimental analysis [5, 19], 

350 MeV < g 0 < 600 MeV . (17) 

Since for the 2 —y 2 processes all momenta of the final particles are fixed by the momenta of 
the initial particles, the above constraint on photon frequency merely defines the region in 
the (Mm, Mm) space where the experimental limit [5, 19] is applicable. 

The cross section for the process + n —> $ + 7 reads 

1 - 2 

= 647rM„ 2 M|u 2 l ' M| 

3 One more requirement on the quantity called asymmetry to be discussed in Sec. 4 is automatically 
fulfilled. 
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where 


2Mq, M n 


tn — U = V 


Finally, we obtain 


n— > 5 > 7 


(My + M n y 
e 2 y 2 P 2 F|(0) 


((M n + M*) 2 - Ml) . 


32 7 T v AMI M x A * 


1 - 


Ml 


(M n + AAp) 


Similarly, the cross section of <f> + n —> T + 7 looks as 

1 eV/3 2 F 2 (0) 


n—>-\I/ 7 


32 7 r v 4 Ml M\ A 4 


1 - 




(Af n + M$) 


1 + 


Mn 

Mu 


The present lower limit on the lifetime of the neutron-decay mode in question is [19, 20] 

T n ^v > 2.8 x 10 31 yr , (18) 


which is applicable in our case while the dark matter masses obey the constraint (17). 
Applying eq. (14), in Fig. 1 we show contours of the constant neutron lifetime of a neutron 
(thin lines) with respect to induced neutron-decay processes < h(T) + n —> T(<h) + 7 . They 
have been calculated for the realistic set of parameters A = M x = 1 TeV and y = 1 without 
any cuts on the phase space. As we explained above, the present experimental limit on this 
process can be applied only within regions shown in violet (light grey) color. In this case, 
one can obtain the current limit on the characteristic scale of the process (A 2 Mx/y) 1 ^] the 
corresponding bounds are shown in these regions by thick lines. Outside shaded blue (dark 
grey) and violet (light grey) regions on this and the subsequent similar plots, the stability 
requirement (4) is not satisfied. 

Note in passing, that applying LHC bounds obtained in [11] is not quite straightforward 
because the couplings ( a which enter (11) are not limited directly from these searches. Thus, 
smaller values of A and M x may be allowed. However, in this paper we will use A = M x = 
1 TeV and y = 1 as a reference set of parameters for numerical estimates. 


4 Scattering processes ^(4>) + n —> 4>(il7) + e + e 


This is a 2 —)• 3 process induced by couplings (13), and (15) through the exchange of a virtual 
neutron and with emission of virtual photon producing an electron-positron pair. With p + 
and p- being the 4-momenta of an outcoming positron and electron, the matrix element is 


ie 2 y/3 F 2 (q 2 )1 + 75 p - M n 


A> c (p 


P 2 - Ml 


(T^q u n (p n ) 47 (p + ) 7 7/> (P-) • 


2 M n A 2 M x q 2 



, GeV 



Figure 1: Contours (thin lines) of constant lifetime (in years) of a neutron with respect to 
the processes 'Fn —>■ $7 and ( I y n —s- T 7 , assuming equal number densities of the two dark 
matter components and parameters A = M\ = 1 TeV and y = 1. Present experimental 
bounds are applicable in the violet (light grey) regions on the plot. Thick lines in these 
regions show the limits on the quantity (A 2 Mx/y) ' in GeV. 

Here q = p + + p~, p = p n — q and T(p^), n(p n ), ijj(p + ), V’(P-) are wave functions of T, n, 
e + , and e~, respectively. 

For the squared matrix element of the 2 —)• 3 process averaged over spins of two incoming 
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fermions we obtain 


(q 2 ■ p n p + ■ qpq, + q 2 • p+pq, ■ qp r 


e*y 2 (d 2 F 2 {q 2 ) 4 1 

1 1 4 Ml M|AV(g 2 -W 

- 4 qpq, • p n p+ • p n p+ + Ap n p + ■ p n p + ■ p n p x] , + 2p n p + ■ qp n ■ qpq. 


+ q 2 M 2 qp xV - q 2 M 2 p n pq, - 4p n p + ■ qp n • p n pq, - 2p + pq, ■ qp n ■ qp n 
+ Aqp n ■ P + P<y ■ p n p + - 2q 2 ■ P + P<y ■ p n p + + 2p n p x y ■ qp n ■ qp n - q 2 • qp n ■ p n p xV ) . 


( 19 ) 


In what follows, it is convenient to describe the final state in terms of energies of the outgoing 
visible particles by choosing, say, positron energy E + and the sum of positron and electron 
energies E. Then for the scattering 4/ n —> <he + e _ , to the leading order in dark matter 
particle velocity v <C 1, one should make the following substitution in eq. (19) (in both the 
center-of-mass and the laboratory frames) 


PnPv = M n Mm , p n p + = M n E + , qp n = EM n , 
p + p<y = E + M x y , qpm = EM x y , qp + = EE + , 

q 2 = 2EM - M 2 + M'l , q 2 - 2qp n = 2EMq, - M 2 + M|, 


where we introduced the notation M = M n + Mq,. Finally, we arrive at 

—^ e 4 y 2 f3 2 F 2 (q 2 ) M n Mq, 

1 hr- e Ml A 4 (,2 _ 2qp n f 

x [q 2 (2 E + (E - E+) - M 2 ) + 2 M 2 (. E\ + (E - E + ) 2 )] . 


The expression for the differential cross section looks as follows [20] 

1 


da = 


4/(27t) 3 16s 


\M\ dm 2 2 dm 


2 

23 ) 


where I = \J ( p n p ^) 2 — M%Mq, ~ M n Mq,v is a flux factor, and in the nonrelativistic limit, 
one has a/s = M. The invariant masses of outgoing pairs (let the subscripts “1” and “2” 
refer to the visible particles and “3” to the dark matter) in the nonrelativistic limit get 
reduced to 


= M 2 + rn 2 — 2MEi , m\ 2 = 2 ME — M 2 + 

2MdE2MdE\. 


dm\ 2 dm^ = 


The energy E is conhned within the interval 

(mi + m 2 ) 2 + M 2 — ml 
2M 


< E < M — m 3 


( 20 ) 


( 21 ) 
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and E\ is within the interval 


where 


and 


ml + M 2 — ( 


m 


2 \max 
23/ 


2 M 


< Ei < 


M 2 + mf 


m 


23 J min 


2 M 



max 

min 


2E* 2 El - m\ - ml ±2 




7-!*2 2 

-^3 — m 3 

p* 2 2 

n / 2 

o 77* 77* 

Z_C / 2 -^3 


M 2 ((M - £) 2 - m 2 ) 

2 ME - M 2 + mf ’ 

(2ME - M 2 + m 2 + m 2 - m 2 ) 2 _ 2 

A{2ME - M 2 + m 2 ) m2 ’ 

(—M 2 + 2ME + ml + mj - ml)(M 2 - ME - mf) 
2MEEM\E~-Vrr? 3 


( 22 ) 


(23) 

(24) 

(25) 


For the process under discussion, let subscript “1” refer to the positron, and replacing in 
the above formulas E\ with E + , we obtain the differential cross section 


da 


1 

1287 


\M\ 2 dEdE + 


which must be integrated over the region defined by eqs. (20)-(25). 

For the process &(p<s>)n(p n ) T(p^) e + (p + )e~(p_), one has the same expression (19) 

multiplied by a factor of 2 due to one less number of initial fermions and makes the replace¬ 
ment 


p n p + = M n E + , qp n = M n E , p n pq, = M n (M - E) , 

<? 2 = _ M 2 + M 2 + 2ME 5 p+pxv = 1 M 2 _ 1 M 2 _ 5 

= M 2 — — ME , g 2 — 2 qp n = —M 2 + + 2 E(M — M n ) , 

where M = + M n . 

The current best limit [20, 19] for neutron decay in the mode n —y ue + e~ is: 

T n ^ ve+e - > 2.57 x 10 32 yr . 

It has been obtained from the analysis of experimental data with imposing the following cut 
on the total energy of leptons [5, 19] 

500 MeV < E < 850 MeV (26) 
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and assuming that the asymmetry is small, 

A <0.5. (27) 

The latter quantity characterizes the directional asymmetry of energy release in the Cherenkov 
detector. The asymmetry is maximal, A — 1. for collinear particles and equals zero for such 
a decay, where the particles go in opposite directions. Let us stress, that this quantity counts 
not all the particles, but only those which release the energy inside the Cherenkov detector, 
and accounts for them with weights proportional to the energy release into the Cherenkov 
radiation. 

In our case of the electron-positron pair the weights are identical. For the decay n —> 
v e + e~ , all 3-momenta of the outgoing particles are in a decay plane. All three particles 
are relativistic, so the Cherenkov angles for the electron and positron are identical and the 
energy conservation gives for the sum of the particle energies 

E u + E-1 _ + E_ = M , (28) 

where M is the neutron mass. Then the asymmetry defined in [5, 19] is just 

A = + n+n ") ’ ( 29 ) 

where n-t are unit 3-vectors along the direction of the outgoing positron and electron, re¬ 
spectively. Introducing the reference axis along the 3-momentum of the neutrino, one defines 
corresponding transverse and longitudinal parts of the electron and positron momenta. Ob¬ 
viously, the transverse parts of electron and positron momenta are equal in magnitude but 
of opposite directions 

P+ = -P-, (30) 

while the longitudinal parts (momentum projection on the chosen axis) sum to zero, 

pI + pI+P 11 - =0. (31) 

For the relativistic electron and positron, one has 

E'l=pi 2 + pi 2 (32) 

and for relativistic neutrino with the chosen axis pi > 0 and 

pi = E u . (33) 
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Then, the asymmetry (29) reads 


, . - , , , P+Pl , Kit 

^ ' TT' TT' ' TT' TT' 

Fj -l JOj— F/+ Fj— 


The differential decay rate is given by 

dT = 1 1 


\M\ dm 2 l2 drn 2 2:i 


(2 7r) 3 32 M 3 

Introducing the sum of the electron and positron energies 

E = E+ + E_ 


(34) 


(35) 


one obtains for the phase space measure (20) (where E\ stands for E + ) that ranges (21) and 
(22) are reduced to 

— <E <M E- — <E + < — . (36) 

2 2 2 

Two independent variables, e.g. E and E +1 fix all the others, which can be found by solving 
Eqs. (28), (30), (31), (32) under condition (33). The results read 


= E - E + , 

E V = M -E, 

E(M - E+) - M 2 /2 


P+ = 
P- = 


M-E 

E(M -E) + EE+ - M 2 /2 
E — M 


_ M ( E ~ M/2)(2E + - M)(E -E + - M/2 ) 
P ~ ~ P+ ~ (E - M) 2 ' 

Putting the solutions above into (34), one obtains for the asymmetry 


(37) 

(38) 

(39) 

(40) 

(41) 

E + (E -E + ) )■ (42) 

The cut adopted in [5, 19] A < 0.5 implies a positive value of the second term in parentheses 
in Eq. (42). It slightly increases the lower limit for E and, thus, reduces a little the triangle 
integration region in (36). 

To adopt the same cuts on asymmetry A in the case of the 2 —> 3 process Tn —» 4>e + e _ 
one can treat it in the nonrelativistic regime as a decay of the particle of effective mass 

M « M n + M* . 

Then the following formulas from the previous considerations must be modified as follows: 


A= - 
2 


E + (E + -E) + M(E - M/2 ) \ 
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• Instead of the massless neutrino, the outcoming dark matter particle $ is massive, so 
its 3-momentum (we use the same notations) instead of (33) obeys 


P\ 


2 + Ml = El 


(43) 


The region of integration in Eq. (36) 


M 2 - Ml 


2 M 


< E < M - M<j>, 


-(E - yJ{M - EY - MJ) < E + < -(E + yJ{M-E)*-Ml), 
Longitudinal momenta are 

E(M - E+) - M 2 f 2 + Ml /2 


P+ = 


P- = 


V(M - Ef - M 2 
E(E + + M — E) — M 2 /2 + Ml/2 


y/(M - E) 2 - M, 


and the transverse momenta read 


j.2 _ ±2 _ [M(E+ -E + M/2 ) - Ml/2][(E - M/2)(M - 2E+) + Mf/2] 
P - ~ P+ ~ (. M-E) 2 -Ml 

The asymmetry (42) must be replaced with 

E + (E + - E) + M(E - M/2) + M|/2 


A = - 1- 
2 


E + (E — E/ 


(44) 


(45) 

(46) 


(47) 


(48) 


Similar formulas with evident replacements and M —> M„ + M$ are applicable 

for the description of the twin process <3 )n —> Te + e _ . 

For the original process n —> ue + e ~, assuming the momenta-independent matrix element, 
the cuts (26) and (27) select a 0.3278/0.3904 part of the phase space. In Fig. 2, we show 
contours of the constant lifetime of a neutron (thin lines) with respect to induced neutron 
decay processes <h('L) + n —> T(4>) + e + e“. They have been calculated without any cuts for 
A = Mx = 1 TeV and y — 1. The current limits on these processes can be applied only within 
the region shown in violet (light grey) color. They are distinguished by the corresponding 
kinematics of the process and applied cuts (26), (27). In this case one can obtain the current 
limit on the characteristic scale of the process A; the corresponding bounds are shown in 
these regions by thick lines. 
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, GeV 


2.8 



Figure 2 : Contours (thin lines) of constant lifetime (in years) of a neutron with respect to 
the process 'h($)?i —* < f>('h)e + e - ; we set A = Mx = 1 TeV and y — 1. Present experimental 
bounds are applicable in the violet (light grey) region on the plot. Thick lines in this region 
show the limits on the quantity (A 2 Mx/y) 1 ^■ 

Now let us consider the asymmetric case when number densities of T and $ are different, 
77 7 ^ 1 , see eq. (3). As an example, below we consider opposite cases of asymmetry: 77 = 0.01 
and 77 = 100, which correspond to $ or $ dominance, respectively. Note that in this case 
the allowed mass intervals are different from that of the symmetric case. Namely, mass 
of the dominant component is fixed in the very narrow region around 5M p /2, while the 
subdominant component can have mass which is determined by the condition (4). In Fig. 3, 
we show expected lifetimes of a neutron with respect to the processes $(T )+?7 —s- 'k(4>)+e + e _ 
for the cases of and T dominance calculated for the same set of parameters as we described 
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Figure 3: Lifetime of a neutron with respect to IND n — y e + e for 77 = 100 and 77 = 0.01; we 
set A = M x = 1 TeV and y — 1. The current limits are applied in shaded regions. 

previously. The current limits on this process are applicable in the shaded regions on these 
figures and they (almost uniformly over these regions) result in A > 35 GeV (77 = 100) and 
A > 40 GeV (77 = 0.01) for region 1 and A > 86 GeV (77 = 100) and A > 75 GeV (77 = 0.01) 
for region 2 . 

Similar plots for the processes $(T) +77 —y T($) + 7 are shown in Fig. 4. Here one can 
obtain the following limits on A: for 77 = 100, we have A > 215 — 228 GeV, and for 77 = 0.01, 
we obtain A > 185 — 200 GeV depending on the mass of the subdominant component. 
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Figure 4: Lifetime of a neutron with respect to the processes $ + n —> 4/ + 7 (for 77 = 100) 
and 'L + n —* <f> + 7 (for y = 0.01); we set A = M x = 1 TeV and y — 1. 


5 Processes + /V —> + 2 mesons 


Within the chiral perturbation theory, the IND processes with two mesons in the final state 
arise in the 1 /f 2 order due to the following terms in the low-energy effective lagrangian 


P Cl 
2/2 


C 


17T 





(t[ + 'K + 



y/3ri7T 0 + 



£271 = ~Yj2 {y§' K v + K ° K ) $ ^ c pr 

^(tt 0 ) 2 + 2 K°K° + K + Kj + h.c., 


(49) 


(50) 


17 






with parameter c\ related to the model parameters as follows from matching eqs. (11) and 
(12) to eqs. (68) and (69) 

= y 

Cl M x A 2 ' 

Details of the derivation are presented in the Appendix A for completeness. Below, we work 
in the limit of exact isotopic invariance, neglecting the proton-neutron and charged-neutral 
pion mass differences, 

M n = M p = M x , m n + = m T o = m n , '%+ = rri^o = mx ■ 

Two types of diagrams contribute the processes: one of them follows from lagrangian (50) 
and the other comes from one-meson lagrangian (49), while the second meson is radiated 
from the nucleon leg; see eq. (72). 

For the dud operator, we have the following possibilities for induced decays, which we 
classify here according to the number of tree-level Feynman diagrams contributing the cor¬ 
responding process: 

• one-diagram processes p —> K°K + , n —> K°K°, and n —> K + K~; 

• two-diagram processes p —> 7i°tt + , n —> n + n~-, 

• three-diagram processes n — > pn 0 , p —$■ r]it + , n —> pp, and n —>• 7r°7r°. 

One-diagram processes. The Feynman diagram for the process T + p —> <F + K°K + 
is presented in Fig. 5 Averaged over spins of the two initial fermions, the squared matrix 



Figure 5: The Feynman diagrams for the process T + p —> <F + K°K + . 


element of this process reads 


\M \ 2 


pc\ 

“g Jl PnP vt, 


(51) 
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and by a factor of two bigger for <f> + p —> T + K°K + ; hereafter, p n refers to the 4- 
momentum of the nucleon participating in the corresponding process. For the first process, 
in the laboratory frame one has to the leading order in dark matter velocity 

PnP^> = Mn Mm (52) 

and when integrating over the phase space adopts the formulas (20)-(25) with 

M = M n + Mm , mi = m 2 = rrix , m 3 = . (53) 

Instead, for the second process we have 

M = M n + Mm , PnP<s> = M n (M — E) , mi = m 2 = m K , m 3 = Mm- (54) 

Averaged over spins of the initial two fermions, squared matrix element of the process 
+ n —> $ + K~K + reads as (51) and by a factor of two bigger for $ + n —> + K~K + . 
Further, in the laboratory frame one can use eqs. (52), (53) and eq. (54) for the first and 
second processes, respectively. The same sets of formulas work for the processes + n — > 

<F + K°K° and $ + n —> + K°K°, respectively. 


Two-diagram processes. To describe this class of processes it is convenient to introduce 
the following notations 


I(p 1 ,p 2 ) = 2pip 2 -pj, (55) 

J(Pi,P 2 ,P 3 ) = 2pip 3 -p 2 p 3 -p\ -pip 2 , (56) 

K (pi,p 2 , p 3 , Pa) = PiPz ■ P2P a + PiPa • P2P3 ~ P1P2 ■ P3P4 ■ (57) 


The Feynman diagrams for the process T + p —> + 7T + 7r° are presented in Fig. 6 . The 

squared matrix element of this process, averaged over spins of initial particles is 


\M \ 2 


(D + F) 2 cl/3 2 M 2 N 

/ 4 

x / J{PN,P*,P*+) J(p N ,P4>,Pno) 
V I 2 (p N ,Pn +) I 2 ip N , Pit 0 ) 


2K (p N , pm , p n +, PirQ ) \ 

I(PN,Pn+)I(pN,Pw°)J ’ 


(58) 


where D = 0.8 and F = 0.47 (see the Appendix A). In the laboratory frame, one has 


p N pm = M N Mm , p N Pn° = M N (E - E x ), p N Pw+ = M N E X , 

p n opm — (E — Ex) Mm , p-n+p^ = E x Mm , p-^+Pn 0 = EM - m 2 n + ^ (. M| - M 2 ) , 
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Figure 6: The Feynman diagrams for the process T + p —> + 7T + 7r°. 

and adopts eqs.( 20 )-( 25 ) with 

M = Mm + m i = m2 = m n , m3 = . 

For $ + p —> T + 7r + 7r°, one obtains for the squared averaged matrix element ( 58 ) but a 
factor of two bigger. In the laboratory frame, one finds 

p N pm = M n (M - E) , p n opm = ^ (M 2 - Af|) - E X M , 

PNPnO = M n (E -El), Pn+p-rrO = EM - m 2 + ^ (Ml - M 2 ) , 

PnPtt+ = M N E l , p n +pm = ^ (M 2 - Ml) - (E - E x ) M , 


M = Mm + M n , rri\ — rn 2 = , m 3 = . 

The predictions of the proton decay with 7r + 7r° final state are presented in Fig. 7 as contours 
of the constant lifetime for the symmetric case rj — 1 . Again, here and below we fix A = 
Mx = 1 TeV and y — 1 and impose no cuts in the phase space. 

Another two-diagram IND process is T + n —> <f> + 7r + n~. Corresponding Feynman 
diagrams are presented in Fig. 8. The squared matrix element of T + n —y + 7r + 7r~ 
averaged over spins of the initial particles takes the form 

1 ( A d\2 ORM m M NP^P^~ | od2 m 2 J(PN,PV,Ptt-) 

-{A - B) p N pm + 2B{A - B) + 2B M N — 

" \PNiPit ) 


where 


(F + D)f3ci 
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Figure 7: Contours of constant lifetime (in years) of the nucleon in the symmetric case with 
respect to IND process with 7r + 7r° in the final state; we set A = Mx = 1 TeV and y — 1. 

In the laboratory frame, one has the same expression as (64), (65) for pions. For $ + n —)• 
T + 7 r + 7 r _ , one has a factor of two bigger squared averaged matrix element and the same 
expressions in the laboratory frame as (66) and (67) for pions. 

Three-diagram processes. The Feynman diagrams for the process T +p —» $ + tt + // are 
presented in Fig. 9. The squared matrix elements of the processes T + p —» $ + 7r + ?/ and 
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Figure 9: The Feynman diagrams for the process T + p —> <f> + ir + r] 


T + n —y <f> + 7r°77, averaged over spins of the initial particles have the form 


2 o/i/f2 ( ry'2 J{PNi P'ih Prr+) J{PN i P^ft Pi~i) o ^ {PN t P^f t P-K +1 Pr]) 


I A A\ z — O I\j2 I u 2 V-^ JV ’ rv > / , r ,i - \j-jv j-?// , 0 uri - *■ vr.iv > r v ■ irrp ( 

m H / 2 (P«,P, + ) U(P«,P,) + /(p»,p,*)/(p»,p,)J 

+ i(A - S - C) 2 PkP „ + 2 B(A - B - op^i + 2C(A - B - C)^i 
2 I{PN,Pn+) I{PN,P V ) 


where 


4 _ y/6Pa D _ [3(D + F)j3ci ^ _ (3 F- D)f3ci 

2/2 ’ V 2 /2 ’ v ^/2 


for T + p —> <f> + Ti + r) and 


_ ^3/9c, _ V3(D + F)/3ci (3F - D)f3ci 

2/2 ’ 2/2 ’ 2y / 3/ 2 
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for \l/ + n —> $ + TT 0, f]- In the laboratory frame, one has 

PnP* = M n M^ , pypr, = My(E - Ei) , PnPtt+ = M N E 1 , 

PnP v = (E- Ex) M n , 2p n +p v = M ( 2 E - M) + M% - - m\ , p n +py = E { M xV , 

and utilizes eqs. (20)-(25) with 

M = My + M n mi = , m 2 = m n , m 3 = . 

For processes <P + p —> T + n + p and <h + n —> T + 7r°?7, one multiplies the above 
expression for the squared averaged matrix element by a factor of 2 and makes the following 
substitutions 


Zp-n+pm = M (M + 2£’i - 2 E) - , 

/Jv/Jvj, = M n (M - E ), 2p, v p r] = M (M - 2F 1 !) - M| + - m 2 , 

P7VPTT+ = M N Ex , p N p r , = (E - E x ) M n , 2 p n +p v = M (2 E - M) + Af| - , 

and uses eqs. (20)-(25) with 


M = M<i, + , mi = , m 2 = nr,,, m 3 = . 


Predictions for the proton lifetime for the 7r + ?7 final state and neutron lifetime for 7i 0 r] final 
state are presented for the symmetric case in Fig, 10 and Fig. 11, respectively. 

The squared matrix element for the processes + n —> $ + 2 77 and + n —> $ + 27r° 
averaged over spins of the initial particles have the form 


\M\ 2 =2 (A - B) 2 p N p^ + AB{A - B ) ^ PlP \ + AB(A - l>) U " v P ' lP " v 


2 B 2 M 2 n 


J(Pn,Pv,Pi) 


+ 


I(Pn,Pi) 
J(PN,P9,P2) 


I 2 (Pn,Pi) I 2 {Pn,P2) 
where p\ and p 2 are momenta of outgoing mesons and 


+ 2 


I(Pn,P-2) 
K(pn,P*,Pi,P2) 
I(Pn,Pi)I(pn,P2) 


_ 3f3ci _ (3F - D)f3c 1 

“ 4/2 ’ “ 2/ 2 

for T + n —> <F + 2r/ and 

Bex (D + F)(3ci 

4 / 2 ’ 2/2 

for if + n —t + 2jr°. In the laboratory frame one has 


(61) 


(62) 

(63) 


PnPi — M N Ei , pvp 2 — My(E - , 2p,p 2 — M (2 E - M) + A !,) - 2 m) , 

PnP2 - ( b’ - /-.'i j M n , - A f v A /a., p,pq, - hi A /,,, 


(64) 
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Figure 10: Contours of the constant lifetime (in years) of a nucleon in the symmetric case 
with respect to IND with tt + t] in the final state; we set A = Mx = 1 TeV and y — 1. 

and adopts eqs. (20)-(25) with 

M = My + , m 3 = , and m i = m 2 = . (65) 

For the averaged squared matrix elements of $ + n —> 'F + 2rj and <f> + n —> 'h + 27 t°, one has 
the same expression (61) multiplied by a factor of two. In the laboratory frame, one finds 

2pip^ — M (M + 2Ei - 2 E) - , p N p x = M N E X , 

2pq,p2 = M (M — 2Ei) — , p N p\v — Mq,(M — E ), (66) 

PnP2 = (E - Ei) M n , 2pip 2 = M (2 E - M) + - 2, 
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Figure 11: Contours of the constant lifetime (in years) of a nucleon in the symmetric case 
with respect to IND with tt°tj in the final state; we set A = Mx = 1 TeV and y — 1. 

and adopts eqs. (20)-(25) with 

M = + M n , m 3 = Mq, , and rri\ = m 2 = m ffl . (67) 

In Fig. 12 we present the predictions of neutron lifetime for the final state 7 r° 7 r°. 

Finally, to illustrate a dependence of the obtained predictions on the value of nonspecified 
asymmetry between and $ populations, 77 , we present in Figs. 13 and 14 the estimates 
of the nucleon lifetime for two opposite cases of large asymmetry r) = 100 and 7] = 0 . 01 , 
respectively. As one observes, the predictions of nucleon lifetimes within hylogenesis model 
can reach values around 10 32 — 10 33 year which looks quite promising for future experiments 
such as Hyper-Kamiokande [1, 16] or DUNE [17]. 
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M*. GeV 

Figure 12: Contours of the constant lifetime (in years) of a nucleon in the symmetric case 
with respect to IND process with 7r°7r° in the final state; we set A = Mx = 1 TeV and y — 1. 

The obtained predictions for double meson channels are, in general, only by an order of 
magnitude weaker than those for single-meson channels (which can be as low as several units 
of 10 31 yr [9, 10] for the same set of parameters). Note that the double meson signatures 
are predicted for the proton decay in the context of grand unified theories [8] as well as for 
dinucleon decays such as pn —y 7T + 7r°, for instance, in supersymmetric models with R-parity 
violation; see e.g., [7]. Searches for the latter type of processes had been performed by the 
Frejus experiment [4] and recently by the Super-Kamiokande collaboration [14]. The most 
stringent limit for the lifetime of the process pn —> 7r + 7r° per oxygen nucleus is found to be 
Tpn^n+nO > 1.70 x 10 32 year. It has been obtained by making use of the expected kinematics 
of dinucleon decay. In particular, the angular distribution between outgoing pions exhibits 
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Figure 13: Contours of the constant lifetime (in years) of the nucleon in the asymmetric 
case, 77 = 100 , with respect to the two-meson processes: a) K°K + ; b) 7 / 77 ; c) 7 r° 7 /; d) 77 + 77 ; 
e) 7 r° 7 r°; g) 7 T+ 7 T 0 . Numbers for other processes with kaons are similar to a) while for process 
with 7T+7T - are similar to g). We set A = Mx = 1 TeV and y — 1. 

a maximum for events with back-to-back topology and the distribution over momentum of 
7T° has a pronounced peak around nucleon mass. Because of this specific kinematics, the 
Super-Kamiokande result cannot be directly applied to the 1ND process with two pions in 
the final state. However, one can show that for some combinations of masses of dark matter 
particles, 4 the distributions over momenta of outgoing mesons also have maxima at 0.5 — 1 
GeV. This signature can be very helpful in discriminating the IND process from the main 
background which is the double pion production by atmospheric neutrinos. 

The double meson channels provide additional signatures of the hylogenesis model, which 
will help to pin down the relevant model parameters once the signal is found. Indeed, even 

4 In particular, when their mass difference is large. 
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Figure 14: Contours of the constant lifetime (in years) of a nucleon in the asymmetric case, 
r) = 0.01 , with respect to the two-meson processes. Other notations are the same as in 
Fig. 13. 

the masses of dark matter particles cannot be unambiguously extracted from a single-meson 
event, because the initial nucleon momentum is not fixed in a real experiment (the nucleon is 
not at rest); hence, the single mesons are not monochromatic. A joint analysis of single and 
double meson events can help to resolve the parameter values. Generally, one anticipates 
that having more than one observable particle in the final state gives more opportunities for 
background reduction in the future experiments. 

The two-meson channels even can help to discriminate between proton decay and induced 
proton decay, which may be challenging is some situations. In particular, if single pions are 
registered at sub-GeV range (say, below 500 MeV), an observation of nmlti-pion events with 
higher total energies would favor the proton decay over the induced proton decay in a model 
where the kinematics constrains the amount of energy allocated to pion at sub-GeV range. 




6 Conclusions 


Summarizing, in this paper we calculated the cross sections of several IND processes for the 
hylogenesis model of dark matter. They include the processes of mimicking neutron decays 
n — y z/y and n —> e + e~. Applying current best limits on the neutron lifetime with respect to 
the processes n —> z/y and n —> e + e~ and taking into account the kinematics of the processes 
which were used in the experiment, we obtained constraints on the parameter space of the 
model. They are considerably weaker than the bounds obtained using the results of the 
searches for events with a high pt jet and missing energy signature at LHC experiments. 

Also, we calculated cross sections and lifetimes corresponding to IND processes with two 
pions in the final state. Searches for such kinds of signatures have not been performed yet 
and present an interesting possibility to further explore the hylogenesis model. We found 
that with the current bounds from the LHC data, the model allows for a lifetime of IND 
such as p — y 7r + n° or n —> at the level of 2 x 10 32 yr. 

Note in passing, that by the time the new generation of experiments looking for nucleon 
decay will be in operation, more data from Run2 of the LHC allow an improvement of the 
collider sensitivity to hylogenesis with respect to the analysis [11]. 

The work was supported by the RSCF grant 14-12-01430. 


A Couplings to baryons and mesons 


The interaction lagrangian of the type (11) with the three light quarks qi = u, q 2 = d, q 3 = s 
in terms of two-component spinors (the relevant are right-handed parts of the Dirac spinors) 
has the form [10] 


Hint = Tr (C O ) + h.c., Oij 


- $ qk%qiR qiR&R, 


( 68 ) 


where 


C = 


( 

\ 


_£2_ I £3 

V6 ^ V2 

0 

0 


0 


C2 _£ 3 _ 

vT V2 


Cl 


0 \ 
0 

~\/i C2 / 


(69) 


The couplings c t are introduced as couplings to the three-quark states which form the eigen¬ 
states of the strong isospin operator. 
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Using the chiral perturbation theory one can obtain [8] the corresponding interaction 
lagrangian for baryons 

C IND = PTr($£CgB R * R )+h.c. 7 


where £ = exp {iM./f) and 


/ _2_ + rU 
1 Ve ^ V 2 


M = 


7T 


K+ \ 


\ 


7 r 

K~ 


JL _zh ts 0 

s/2 ^ 

K° 


b) 


and baryon fields leaving only a neutron and proton 


Br = 


( 0 0 p R \ 

0 0 n R 


\ 0 0 0 / 

Expanding to linear order in meson fields we find (hereafter, in terms of the Dirac fermions) 


2 j ^ _ 

£| 7 T = y v b 6 ' Cl ^ -y -nfi 'r/ + -j=n R ir” - p R it 

c 2 V3 c 3 \ /c 3 \ - 0 

vr + viJ PRA + . 

Expanding to the second order in 1//, one obtains 

£2, = jL (-V;i • * C PJJ $ + A 32 • f c riH'[') - I'.C. 


( 70 ) 


+ h.c. 


( 71 ) 


where 


t 4 3 i = - ci (y6n 77 + K°K j 


3 Vs \ (Vs 1 , 0 

§ C2 + — C 3 \K--r, + — C2 + -C3 \K -2T 


^ - -|c 3 ] 


t4 3 2 — 


— Ci ^7r + 7r + bj 2 — VSr]7i° + - (7r 0 ) 2 + 2 K°K° + K + K ^ 

Vs 3c 3 ^ + (3 V3 \ - q / Vs c 3 1 (, q 

_ f2 + - U I++ - C2 - T C3 A -c 2 - ? |A .. 


30 



Finally, for completeness let us remind [8] here the interaction lagrangian of baryons with 
mesons to the leading order in derivative expansion, which has the form 


o rp _ 7W _i_ rp 

C = ——-(p7 M 7 5 p + h7 M 7 5 n) d^r] H-(p7 M 7 5 p — d^n 0 


V6f 


V2 


D + F 

H - 7 — (<9 /i 7r + p7 / Vn + <9 M 7r“h7 M 7 5 p) , 


(72) 


/ 


where D = 0.8 and F = 0.47. 
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